Lagrangian Structure Functions in Turbulence: A Quantitative Comparison between 

Experiment and Direct Numerical Simulation 
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A detailed comparison between data from experimental measurements and numerical simulations 
of Lagrangian velocity structure functions in turbulence is presented. By integrating information 
from experiments and numerics, a quantitative understanding of the velocity scaling properties over 
a wide range of time scales and Reynolds numbers is achieved. The local scaling properties of the 
Lagrangian velocity increments for the experimental and numerical data are in good quantitative 
agreement for all time lags. The degree of intermittency changes when measured close to the 
Kolmogorov time scales or at larger time lags. This study resolves apparent disagreements between 
experiment and numerics. 

PACS numbers: 47.27.Eq ,47.27.Gs, 02.50.-r, 47.27.Jv 



I. INTRODUCTION 



the Lagrangian Velocity Structure Functions (LVSF), 
defined as 
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Understanding the statistical properties of a fully 
developed turbulent velocity field from the Lagrangian 
point of view is a challenging theoretical and experimen- 
tal problem. It is a key ingredient for the development 
of stochastic models for turbulent transport in such 
diverse contexts as combustion, pollutant dispersion, 
cloud formation, and industrial mixing ji^ 3 .^ Progress 
has been hindered primarily by the presence of a wide 
range of dynamical timescales, an inherent property 
of fully developed turbulence. Indeed, for a complete 
description of particle statistics, it is necessary to 
follow their paths with very fine spatial and temporal 
resolution, on the order of the Kolmogorov length and 
time scales r\ and r^. Moreover, the trajectories should 
be tracked for long times, order the eddy turnover 
time Tl, requiring access to a vast experimental mea- 
surement region. The ratio of the above timescales 
can be estimated as Tl/t v ~ R\, and the microscale 
Reynolds number R\ ranges from hundreds to thou- 
sands in typical laboratory experiments. Despite these 
difficulties, many experimental and numerical studies 
of Lagrangian turbulence have been reported over the 

,,..„ 5.6.7.8.9 .10.11.12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28 . 

Here, we present a detailed comparison between state- 
of-the-art experimental and numerical studies of high 
Reynolds number Lagrangian turbulence. We focus on 
single particle statistics, with time lags ranging from 
smaller than r n to order Tj,. In particular, we study 



S p (r) = ((5 T vr) = ([v(t + r)-v(t)r), 



(1) 



where v denotes a single velocity component. 

In the past, the corresponding Eulerian quantities, 
i.e. the moments of the spatial velocity increments, have 
attracted significant interest in theory, experiments, and 
numerical studies (for a review see Ref. |35T ) . It is now 
widely accepted that spatial velocity fluctuations are in- 
termittent in the inertial range of scales, for rj <C r <C L, 
L being the largest scale of the flow. By intermittency 
we mean anomalous scaling of the moments of the veloc- 
ity increments, corresponding to a lack of self-similarity 
of their probability density functions (PDFs) at different 
scales. In an attempt to explain Eulerian intermittency, 
many phenomenological theories have been proposed, ei- 
ther based on stochastic cascade models {e.g. multifractal 
description a 36 ' 37 i 38 ), or on closures of the Navier-Stokes 
equations^ Common to all these models is the presence 
of non-trivial physics at the dissipative scale, r ~ rj, in- 
troduced by the complex matching of the wild fluctua- 
tions in the inertial range and the dissipative smoothing 
mechanism at small scalesi 40 ! 41 Numerical and experi- 
29-3 fli 3 ii 32 :P 3 j?:'iervations show that clean scaling behavior for 



the Eulerian structure functions is found only in a range 
10i] < r -C L (see Refl42l for a collection of experimental 
and numerical results). For spatial scales r < 10i], multi- 
scaling properties, typical of the intermediate dissipative 
range, are observed due to the superposition of inertial 
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range and dissipative physics 4^ 

Similar questions can be raised in the Lagrangian 
framework: (i) is there intermittency in Lagrangian 
statistics? (ii) is there a range of time lags where clean 
scaling properties (i.e. power law behavior) can be de- 
tected? (iii) are there signatures of the complex interplay 
between inertial and dissipative effects for small time lags 
T~0(r n )? 

In this paper we shall address the above questions by 
comparing accurate Direct Numerical Simulations (DNS) 
and laboratory experiments. Unlike Eulerian turbulence, 
the study of which has attracted experimental, numer- 
ical and theoretical efforts since the last thirty years, 
Lagrangian studies become available only very recently 
mainly due to the severe difficulty of obtaining accu- 
rate experimental and numerical data at sufficiently high 
Reynolds numbers. Consequently, the understanding of 
Lagrangian statistics is still poor. This explains the ab- 
sence of consensus on the scaling properties of the LVSF. 
In particular, there have been different assessments of the 
scaling behavior 

S p (r) = ((S T vy)~T^, (2) 

mainly due the desire to extract a single number, i.e. the 
scaling exponent £(p), over a range of time lags. 

Measurements using acoustic technique o 10 ' 15 gave the 
first values of the exponents £ (p) , measuring scaling prop- 
erties in the range 10r ri < t < T^. Subsequently, exper- 
iments based on CMOS sensor s 26 ' 28 provided access to 
scaling properties for shorter time lags, 2r r) < r < 6^, 
finding more intermittent values, though compatible with 
Ref. [lCJ. DNS data, obtained at lower Reynolds num- 
ber, allowed simultaneous measurements in both of these 
rangesi 23 ' 29 For I0r n < r < 50t v , scaling exponents were 
found to be slightly less intermittent than those measured 
with the acoustic techniques, though again compatible 
within error bars. On the other hand, DNS dat a 29 ' 34 ' 43 
for small time lags, 2T n < r < 6r^, agree with scaling 
exponents measured in Ref. l26l 

The primary goal of this paper is to critically compare 
state-of-the-art numerical and experimental data in or- 
der to analyze intermittency at both short and long time 
lags. This is a necessary step both to bring Lagrangian 
turbulence up to the same scientific standards as Eulerian 
turbulence and to resolve the conflict between experiment 
and simulations (see also Refs. HHIHIHII) . 

To illustrate some of the difficulties discussed above, 
in Fig. [T] we show a compilation of experimental and 
numerical results for the second-order Lagrangian struc- 
ture function at various Reynolds numbers (see later for 
details). The curves are compensated with the dimen- 
sional prediction given by the classical Kolmogorov di- 
mensional theory in the inertial range 4 ^, Si(j) = Cqct, 
where e is the turbulent kinetic energy dissipation. The 
absence of any extended plateau and the trend with the 
Reynolds number indicate that the inertial range, if any, 
has not developed yet. The same trends have been ob- 
served in other DNS studies 2 ^ and by analyzing the tern- 
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FIG. 1: Log-log plot of the second-order LVSF (averaged over 
the three components) normalized with the dimensional pre- 
diction, i.e. 5 l 2(r)/(eT), at various Reynolds numbers and for 
all data sets. Details can be found in Tables 1 and 2. EXP2 
and EXP4 refer to the same Reynolds number (R\ = 690), 
but with different measurement volumes (larger in EXP4); 
in particular EXP2 and EXP4 better resolve the small and 
large time lag ranges, respectively, and intersect for t/t v ~ 2. 
We indicate with a solid line the resulting data set made of 
data from EXP2 (for t/t v < 2) and EXP4 (for t/t v > 2); 
a good overlap among these data is observed in the range 
2 < t/t v < 8. For all data sets, a extended plateau is absent, 
indicating that the power law regime typical of the inertial 
range has not yet been achieved, even at the highest Reynolds 
number, R\ ~ 815, in experiment. 



poral behavior of signals with a given power-law Fourier 
spectrum^ 6 - 

We stress that assessing the actual scaling behavior of 
the second (and higher) order Lagrangian velocity struc- 
ture functions is crucial for the development of stochastic 
models for Lagrangian particle evolution. Indeed, these 
models are based on the requirement that the second- 
order LVSF scales as S2 (r) cx er. The issues of whether 
the predicted scaling is ever reached and ultimately how 
the LVSF deviate as a function of the Reynolds numbers 
remains to be clarified. 

Moreover, an assessment of the presence of Lagrangian 
intermittency calls for more general questions about phc- 
nomenological modeling. For instance, multifractal mod- 
els derived from Eulerian statistics can be easily trans- 
lated to the Lagrangian framework ; 10 i 23 i 47 ' 48 with some 
degree of successi 10 ' 13 ' 18 

The material is organized as follows. In Section |TT1 we 
describe the properties of the experimental setup and the 
direct numerical simulations, detailing the limitations in 
both sets of data. A comparison of Lagrangian veloc- 
ity structure functions is considered in Section [TTT1 Sec- 
tion IIII Al presents a detailed scale-by-scale discussion of 
the local scaling exponents, which is the central result of 
the paper. Section llVI draws conclusions and offers per- 
spectives for the future study of Lagrangian turbulence. 
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II. EXPERIMENTS AND NUMERICAL 
SIMULATIONS 

Before describing the experimental setup and the DNS 
we shall briefly list the possible sources of uncertainties 
in both experimental and DNS data. In general this is 
not an easy task. First, it is important to discern the de- 
terministic from the statistical sources of errors. Second, 
we must be able to assess the quantitative importance of 
both types of uncertainties on different observables. 

Deterministic uncertainties. For simplicity, we report 
in this work the data averaged over all three components 
of the velocity for both the experiments and the DNS. 
Since neither flows in the experiments nor the DNS are 
perfectly isotropic, a part of the uncertainty in the re- 
ported data comes from the anisotropy. In the experi- 
ments the anisotropy reflects the generation of the flow 
and the geometry of the experimental apparatus. The 
anisotropy in DNS is introduced by the finite volume 
and by the choice of the forcing mechanism. In gen- 
eral, the DNS data are quite close to statistical isotropy, 
and anisotropy effects are appreciable primarily at large 
scales. This is also true for the data from the experiment, 
especially at the higher Reynolds numbers. An important 
limitation of the experimental data is that the particle 
trajectories have finite length due both to finite mea- 
surement volumes and to the tracking algorithm, which 
primarily affect the data for large time lags. It needs to 
be stressed, however, that in the present experimental 
set up due to the fact that the flow is not driven by bulk 
forces, but by viscous and inertial forces at the blades, 
the observation volume would anyhow be limited by the 
mean velocity and the time it takes for a fluid particle 
to return to the driving blades. At the blades the tur- 
bulence is strongly influenced by the driving mechanism. 
Therefore, in the experiments reported here the observa- 
tion volume was selected to be sufficiently far away from 
the blades to minimize anisotropy. For short time lags, 
the greatest experimental difficulties come from the finite 
spatial resolution of the camera and the optics, the im- 
age acquisition rate, data filtering and post-processing, a 
step necessary to reduce noise. For DNS, typical sources 
of uncertainty at small time lags are due to the inter- 
polation of the Eulerian velocity field to obtain the par- 
ticle position, the integration scheme used to calculate 
trajectories from the Eulerian data, and the numerical 
precision of floating point arithmetic. 

The statistical uncertainties for both the experimen- 
tal and DNS data arise primarily from the finite number 
of particle trajectories and-especially for DNS-from the 
time duration of the simulations. We note that this prob- 
lem is also reflected in a residual, large-scale anisotropy 
induced by the non-perfect averaging of the forcing fluc- 
tuations in the few eddy turnover times simulated. The 
number of independent flow realizations can also con- 
tribute to the statistical convergence of the data. While 
it is common to obtain experimental measurements sepa- 
rated by many eddy turnover times, typical DNS results 



contain data from at most a few statistically independent 
realizations. 

We stress that, particularly for Lagrangian turbulence, 
only an in-depth comparison of experimental and numer- 
ical data will allow the quantitative assessment of uncer- 
tainties. For instance, as we shall see below, DNS data 
can be used to investigate some of the geometrical and 
statistical effects induced by the experimental apparatus 
and measurement technique. This enables us to quantify 
the importance of some of the above mentioned sources 
of uncertainty directly. DNS data are, however, limited 
to smaller Reynolds number than experiment; therefore 
only data from experiments can help to better quantify 
Reynolds number effects. 



A. Experiments 

The most comprehensive experimental data of La- 
grangian statistics are obtained by optically tracking pas- 
sive tracer particles seeded in the fluid. Images of the 
tracer particles are analyzed to determine their motion in 
the turbulent flowi&i^ Due to the rapid decrease of the 
Kolmogorov scale with Reynolds number in typical lab- 
oratory flows, previous experimental measurements were 
often limited to small Reynolds numbersi 6 -^ The Kol- 
mogorov time scale at R\ ~ 0(1 3 ) in a laboratory water 
flow was so far resolved only by using four high speed sil- 
icon strip detectors originally developed for high-energy 
physics experimentsi&ii The one-dimensional nature of 
the silicon strip detector, however, restricted the three 
dimensional tracking to a single particle at a time, limit- 
ing severely the rate of data collection. Recent advances 
in electronics technology now allow simultaneous three 
dimensional measurements of O(10 2 ) particles at a time, 
by using three cameras with two-dimensional CMOS sen- 
sors. High-resolution Lagrangian velocity statistics at 
Reynolds numbers comparable to those measured using 
silicon strip detectors are therefore becoming available, 26 

Lagrangian statistics can also be measured acousti- 
cally. The acoustic technique measures the Doppler fre- 
quency shift of ultrasound reflected from particles in the 
flow, which is directly proportional to their velocity 10 ^ 5 
The size of the particles needed for signal strength in the 
acoustic measurements can be significantly larger than 
the Kolmogorov scale of the flow. Consequently, the par- 
ticles do not follow the motion of fluid particles^ and 
this makes the interpretation of the experimental data 
more difficult^ 

The experimental data here presented are discussed in 
much detail in Refs. [26,28. In the following, we only 
briefly recall the main aspects of the experimental tech- 
nique and data sets, whose parameters are summarized 
in Table [J 

Turbulence was generated in a swirling water flow be- 
tween counter-rotating baffled disks in a cylindrical con- 
tainer. The flow was seeded with polystyrene particles 
of size d p = 25/im and density p p = 1.06g/cm 3 that 
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EXP 2 
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0.42 
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0.16 


24 


0.3x0.3x0.3 
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480 


9.6xl0 5 


EXP 3 


815 


0.59 


3.0 


23 


0.54 


0.11 


15 


0.3x0.3x0.3 


80 


500 


1.7xl0 6 


EXP4 


690 


0.42 


1.2 


30 


0.90 


0.16 


24 


0.7x0.7x0.7 


200 


1200 


6.0xl0 6 



TABLE I: Parameters of the experiments. Column three gives the value of the root-mean-square velocity fluctuations v' rms , 
averaged over the three components. The integral length scale L = v' rms 3 /e = 7cm was determined to be independent of 
Reynolds number. Tl = L/v' rms is the eddy turnover time. Nf is the temporal resolution of the measurement, in units of 
frames per t v . The measurement volume was nearly a cube in the center of the tank and its linear dimensions are given in 
units of the integral length scale L. Ax is the spatial discretization of the recording system. The spatial uncertainty of the 
position measurements is roughly 0.1 Ax. Nr is the number of independent realizations recorded (see text). N tr is the number 
of Lagrangian trajectories measured. 



follow the flow faithfully for R\ up to lO 3 ^ The parti- 
cles were illuminated by high-power Nd:YAG lasers, and 
three cameras at different viewing angles were used to 
record the motion of the tracer particles in the center 
of the apparatus. Images were processed to find particle 
positions in three-dimensional physical space; the parti- 
cles were then first tracked using a predictive algorithm 
to obtain the Lagrangian trajectories^ Due to fluctu- 
ations in laser intensity, the uneven sensitivity of the 
physical pixels in the camera sensor array, plus electronic 
and thermal noise, images of particles sometimes fluctu- 
ate and appear to blink. When the image intensity of a 
particle was too low, the tracking algorithm lost that par- 
ticle. Consequently, the trajectory of that particle was 
terminated. When the image intensity is high again, the 
algorithm started a new trajectory. The raw trajectories 
therefore contained many short segments that in reality 
belonged to the same trajectory. It is, however, possible 
to connect these segments by applying a predictive al- 
gorithm in the six-dimensional space of coordinates and 
velocities^ The trajectories discussed in this paper were 
obtained with the latter method which allows for much 
longer tracks. 

The Lagrangian velocities were calculated by smooth- 
ing the measured positions and subsequently differentiat- 
ing. A Gaussian filter has been used to smooth the data. 
Smoothing and differentiation can be combined into one 
convolution operation by integration by parts; the con- 
volution kernel is simply the derivative of the Gaussian 
smoothing filter The width of the Gaussian kernel was 
chosen to remove the noise in position measurements, 
but not to suppress the fluctuations, whose characteristic 
time scale is 0(t v ) or above. The velocity statistics have 
been found to be insensitive to the width a of the Gaus- 
sian filter, provided it is between t v /6 and (see also 
below). The temporal resolution of the camera system 
in the experiments reported here was sufficiently high to 
ensure that the fluctuations with time scale greater than 
were well resolved. 

The uncertainty in position measurement, or the spa- 
tial resolution, is directly proportional to the size of the 



spatial discretization determined by the optical magnifi- 
cation and by the size of the pixels on the CMOS sensor. 
Larger magnification gives better spatial resolution but 
also a smaller measurement volume. Indeed the num- 
ber of pixels of the camera sensor array is fixed by the 
chip-size and, at higher speeds, by the imaging rate. The 
dynamic range of the cameras is not sufficient to cover the 
entire range of scales of the turbulence at the Reynolds 
numbers of interest. Therefore, two sets of experiments 
with different magnifications have been performed. The 
former set has high spatial resolution and focuses on 
the small scale quantities, though with a relatively small 
measurement volume (EXP 1,2, 3 in Table J]). Then, in 
order to probe longer times and larger scales, the size of 
the measurement volume in the second set of measure- 
ments was chosen to be slightly smaller than the inte- 
gral scale (EXP4 in Table [TJ) . In this data set, however, 
the uncertainty in position was larger and the short-time 
statistics were severely affected. As a result, in order to 
have experimental data covering a wide range of time lags 
(tjj < t < IOOt,,) at a given Reynolds number, one needs 
to merge data from the two different experiments. This 
could be done at Re\ = 690, by using data from the small 
measurement volume (EXP2) up to times t ~ (6 -j- 7)r^, 
and using data from large measurement volume (EXP4) 
at larger times. The procedure is well justified as the two 
data sets match for intermediate time lags. 

One noticeable difference between experiments and nu- 
merical simulations is the number of independent real- 
izations included in the statistics. While it is difficult to 
have many statistically independent DNS results at one 
Reynolds number, the experimental data usually con- 
tained O(10 3 ) records separated by a time interval of 
about 10 2 Tl. Each of these records lasted for (1 2) T^. 
The variation of the velocity fluctuations calculated from 
the statistics of many records is shown in Figure 2(a) As 
it is clear from the figure, the three components do not 
fluctuate about the same value, indicating the presence 
of anisotropy which does not average away even after 
many eddy turnover times. These effects are introduced 
by the flow generation in the apparatus. In the following 
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the uncertainties in the data sets due to anisotropy were 
estimated by the difference between measurements made 
on different components of the velocity field. 



B. Direct Numerical Simulations 

Nowadays state-of-the-art numeric a 19 i 23 ' 51 i 52 best 
suited for Eulerian statistics is able to reach Taylor scale 
Reynolds numbers of the order of R\ ~ 1000 by using up 
to 4096 3 mesh points^ Such extremely high Reynolds 
number DNS is, however, limited by the impossibility 
of integrating the flow for long time durations, due to 
the extremely high computational costs. In Lagrangian 
studies it is necessary to highly resolve the Eulerian veloc- 
ity field to obtain precise out-of-grid interpolation. The 
maximum achievable Reynolds number, on the fastest 
computers, is currently limited to R\ ~ 600 in order to 
accurately calculate the particle positions and to achieve 
sufficiently long integration times4 ' 19 i 23 ' 27 

Typically, such Lagrangian simulations last for a few 
large-scale eddy turnover times, implying some unavoid- 
able remaining anisotropy at large scales, even for nom- 
inally perfectly isotropic forcing. The simulations ana- 
lyzed here were forced by fixing the total energy of the 
first two Fourier-space shells^ 3 -: E{k\) = J2\k\eh 1^(^)1 
and E(k 2 ) = J2\k\ei 2 1^(^)1' wnere h = [0-5 : 1.5] and 
I2 = [1.5 : 2.5] (the \k\ — mode is fixed to zero to 
avoid a mean flow). The three velocity components can 
instantaneously be quite different: when one of the three 
fluctuates, the others must compensate in order to keep 
the total amplitude fixed (see, for instance, Fig. 2(b) for a 
visualization of this effect). However, by averaging over 
many eddy turn over times - when possible, as for the 



lower-resolution DNS shown in the inset of Fig. |2(b)[ 
, the forcing produces a perfectly statistically isotropic 
flow. As the remaining large-scale anisotropy is the main 
source of uncertainty in the DNS results, we will estimate 
confidence intervals from the difference between the three 
components. 

In the simulations, the main systematic error for small 
time lags comes from the interpolation of the Eulerian ve- 
locity fields needed to integrate the equation for particle 
positions, 



X(t)=v(X(t),t) 



(3) 



Of course, high-order interpolation schemes such as third- 
order Taylor series interpolation or cubic splines partially 
remove this problem. Cubic splines give higher interpo- 
lation accuracy, but they are more difficult to use in im- 
plementations that rely on secondary storage i 54 ' 55 It has 
been reported^ that cubic schemes may resolve the most 
intense events better than linear interpolation, especially 
for acceleration statistics; the effect, however, appears to 
be rather small especially as far as velocity is concerned. 

More crucial than the order of the interpolation scheme 
is the resolution of the Eulerian grid in terms of the Kol- 
mogorov length scale. To enlarge the inertial range as 



much as possible, typical Eulerian simulations tend to 
poorly resolve the smallest scale velocity fluctuations by 
choosing a grid spacing Ax larger then the Kolmogorov 
scale r\. Since this strategy may be particularly harmful 
to Lagrangian analysis, here it has been chosen to bet- 
ter resolve the smallest fluctuations by choosing Ax ~ r\ 
and to use the simple and computationally less expensive 
linear interpolation. 

We stress that having well resolved dissipative physics 
for the Eulerian field is also very important for captur- 
ing the formation of rare structures on a scale r »s 77. 
Moreover, as discussed in Ref. [13, such structures, be- 
cause of their filamentary geometry, may influence not 
only viscous but also inertial range physics. 
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FIG. 2: (b) Time evolution of the components of the ve- 
locity fluctuation v'£ (dashed line), v'y (thick black line) and 
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v z ~ (solid line) for EXP2. (a) Time evolution of vf, with 
i — x,y,z, for DNS2. In the inset we show the same time 
evolution for a DNS at a smaller R\ w 75 (obtained with a 
spatial resolution of 128 3 grid points and the same forcing), 
which was integrated for a much longer time. In the latter 
case, the three components fluctuate around the same value, 
showing the recovery of isotropy for long enough time. 
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DNS1 
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1.5 


0.886 
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0.01 


3.14 


2.1 


0.048 


5 


0.012 


512 3 


0.96xl0 6 


DNS2 
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0.005 


3.14 
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TABLE II: Parameters of the numerical simulations. Taylor microscale Reynolds number R\, root-mean-square velocity 
fluctuations v' rms , energy dissipation e, viscosity i>, Kolmogorov length scale r\ = (y 3 /e) 1//4 , integral scale L, large-eddy turnover 
time Tl = L/v' rrns , Kolmogorov time scale t v = {u/e) 1 ^ 2 , total integration time T, grid spacing Ax, resolution JV 3 , and the 
number of Lagrangian tracers N p . 



Another possible source of error comes from the loss 
of accuracy in the integration of Eq. ((3|) for very small 
velocities due to round-off errors. This problem can be 
overcome by adopting higher-order schemes for tempo- 
ral discretization. For extremely high Reynolds numbers 
it may also be necessary to use double precision arith- 
metic, while for moderate R\, single precision, which was 
adopted in the present DNS, is sufficient for accurate re- 
sults (see, e.q.. Ref.il). 

Details of the DNS analyzed here can be found else- 
where^; here, we simply state that the Lagrangian trac- 
ers move according to Eq. ([3]), in a cubic, triply periodic 
domain of side B = 2ir. DNS parameters are summarized 
in Table HJ 



III. COMPARISON OF LAGRANGIAN 
STRUCTURE FUNCTIONS 

Let us now compare the experimental and numeri- 
cal measurements of the Lagrangian velocity structure 
functions directly. Figures 3(a) and 3(b) show a di- 
rect comparison of LVSFs of order p = 2 and p = 4 
for all data sets. The curves are plotted using the di- 
mensional normalization, assuming that £2(1") = Cq€t cx 



,i? A 1 (r/r^) (where we use e 



J L and Tl/t v cx 



R\). Such a rescaling can be generalized as S p (t) cx 

IP D -p/2( T / T ^p/2 4 Both the 2 "d and 4th order 



V rms^X 



ill 



2T mo- 



ments show a fairly good collapse, especially in the range 
of intermediate time lags. However, some dependence 
can be observed both on R\ (see Fig. |3(b)[ ) and on the size 
of the measurement volume (compare EXP2 and EXP4). 
Both effects call for a more quantitative understanding. 



A. Local Scaling Exponents 

A common way to assess how the statistical properties 
change for varying time lags is to look at dimensionless 
quantities such as the generalized flatness 



F2 P (r) 



[5a(r)]f 



(4) 



We speak of intermittency when such a function changes 
its behavior as a function of r: this is equivalent to the 
PDF of the velocity fluctuations S T v, normalized to unit 
variance, changing shape for different 
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FIG. 3: (a) Log-log plot of the second-order structure function 
compensated as R\S2(j)/v' rma 2 vs t/t,, for all data sets, at 
several Reynolds numbers, (b) The same for the fourth-order 
structure function R\Si{r)/v' rms . The solid line is made to 
guide the eye through the two data sets (EXP2 and EXP4) 
obtained at the same Reynolds number in two different mea- 
surement volumes, as explained in Sect. Ill Al 



When the generalized flatness varies with r as a power 
law, F2 P (t) ~ T x ^ 2p \ the scaling laws are intermittent. 
Such behavior is very difficult to assess quantitatively, 
since many decades of scaling are typically needed to 
remove the effects of sub-leading contributions (for in- 
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FIG. 4: Generalized flatness F2 P (t) of order p — 2 and p = 3, 
measured from DNS2, EXP2, and EXP4. Data from EXP 2 
and EXP4 are connected by a continuous line. The Gaussian 
values are given by the two horizontal lines. The curves have 
been averaged over the three velocity components and the 
error bars are computed from the scatter between the three 
different components as a measure of the effect of anisotropy. 
Statistical errors due to the limitation in the statistics are 
evaluated by dividing the whole data sets in sub samples 
and comparing the results. These statistical errors are always 
smaller than those estimated from the residual anisotropy. 



stance, it is known that Eulerian scaling may be strongly 
affected by slowly decaying anisotropic fluctuations^). 

We are interested in quantifying the degree of intermit- 
tency at changing r. In Fig. [4] we plot the generalized 
flatness F2 P (t) for p = 2 and p — 3 for the data sets 
DNS2, EXP2 and EXP4. Numerical and experimental 
results are very close, and clearly show that the inter- 
mittency changes considerably going from low to high r. 

The difficulty in trying to characterize these changes 
quantitatively is that, as shown by Fig. [4j one needs to 
capture variations over many orders of magnitude. For 
this reason, we prefer to look at observables that remain 
0(1) over the entire range of scales and which convey in- 
formation about intermittency without having to fit any 
scaling exponent. With this aim, we measured the loga- 
rithmic derivative (also called local slope or local expo- 
nent) of structure function of order p, S p (t), with respect 
to a reference structure function^ for which we chose the 
second-order ft (r): 



d log (S P (r)) 
dbg (ft (t)) 



(5) 



We stress the importance of taking the derivative with re- 
spect to a given moment: this is a direct way of looking at 
intermittency with no need of ad hoc fitting procedures 
and no request of power law behavior. This procedure,— 
which goes under the name of Extended Self Similar- 
ity^ (ESS), is particularly important when assessing the 
statistical properties at Reynolds numbers not too high 



and/or close to the viscous dissipative range. 

A non-intermittent behavior would corresponds to 
Cp( T ) = p/2< In the range of t for which the exponents 
Cp(t) are different from the dimensional values p/2, struc- 
ture functions are intermittent and correspondingly the 
normalized PDFs o f (S T v / {(Srv) 2 ) 1 ^ 2 ) change shape with 
t. Figures 5(a) and 5(b) show the logarithmic local slopes 



p = 4 
from the 



of the numerical and experimental data sets for several 
Reynolds numbers for p — 4 and p = 6 versus time nor- 
malized to the Kolmogorov scale, t/t v . These are the 
main results of our analysis. 

The first observation is that for both orders p — 4 and 
p = 6, the local slopes Cp( T ) deviate strongly from their 
non-intermittent values £4 — 2 and = 3. There is a 
tendency toward the differentiable non-intermittent limit 
Cp = p/2 only for very small time lags r <C r, ; . 

In the following, we shall discuss in detail the small 
and large time lag behavior. 

Small time lags. For the structure function of order 
(Fig. 5(a)), we observe the strongest deviation 
non-intermittent value in the range of time 
2t^ < t < Gt^. It has previously been proposed that this 
deviation is associated with particle trapping in vortex 
filaments^ This fact has been supported by DNS inves- 
tigations of inertial particles i 14 i 17 ' 29 The agreement be- 
tween the DNS and the experimental data in this range is 
remarkable. For p = 6 (Fig. 5(b) I, the scatter among the 
data is higher due to the fact that, with increasing order 
of the moments, inaccuracies in the data become more 
important. Still, the agreement between DNS and the 
experimental data is excellent. Differently from the p = 4 
case, a dependence of mean quantities on the Reynolds 
number is here detectable, though it lies within the error- 
bars. The experimental data set for p — 6, at the highest 
Reynolds number (R\ = 815), show a detectable trend in 
the local slope toward less intermittent values in the dip 
region, 2 < t/t v < 6. This change may potentially be 
the signature of vortex destabilization at high Reynolds 
number - which would reduce the effect of vortex trap- 
ping. It is more likely, however, that at this very high 
Reynolds number both spatial and temporal resolution 
of the measurement system may not have been sufficient 
to resolve the actual trajectories of intense events^ We 
consider this to be an important open question for future 
studies. 

Larger time lags. For t > (6 -j- 7)t v up to Tl, the 
experimental data obtained in small measurement vol- 
umes (EXP1,2,3), are not resolving the physics, as they 
develop both strong oscillations and a common trend to- 
ward smaller and smaller values for the local slopes for 
increasing r. This may be attributed to finite volume 
corrections (see also Sect. IIII All) . For these reasons, the 
data of EXP1,2,3 are not shown for these time ranges. 
On the other hand, the data from EXP4, obtained from 
a larger measurement volume, allow us to compare exper- 
iment and simulation. Here the local slope of the experi- 
mental data changes slower very much akin to the simu- 
lations. This suggests that in this region high Reynolds 
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FIG. 5: Logarithmic derivatives Cp( t ) or structure functions S p (t) with respect to £2(1") for orders p — 4 (a) and p = 6 (b). 
The curves are averaged over the three velocity components and the error bars are computed from the statistical (anisotropic) 
fluctuations between LVSFs of different components. The horizontal lines are the non-intermittent values for the logarithmic 
local slopes, i.e. £ p = p/2. We stress that the curves for EXP1,2,3 are shown in the time range f < t/t v < 7, while the curves 
for EXP4 (large measurement volume) are shown in the time range 7 < t/t v < 50. 



number turbulence may show a plateau, although the larger uncertainties. The excellent quantitative agree- 
current data can not give a definitive answer to this ques- ment between DNS and the experimental data gives us 
tion. For p = 6, a similar trend is detected, though with high confidence into the local slope behavior as a function 
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of time lag. 

In light of these results, we can finally clarify the re- 
cent apparent discrepancy between measured scaling ex- 
ponents of the LVSFs in experiments^ and DNS,— which 
have lead to some controversy in the literature. 34 ! 43 ' 44 In 
the experimental work>2£, scaling exponents were mea- 
sured by fitting the curves in Fig. [5] in the range 2t v < 
t < 6r^, where the compensated second order veloc- 
ity structure functions reach a maximum, as shown in 
Fig. [U (measuring the fourth and sixth order scaling ex- 
ponents Cp( T ) to be 1.4 ±0.1 and 1.6 ±0.1, respectively). 
On the other hand, in the simulations^ 3 - scaling expo- 
nents were measured in the regions in the range of times 
10t,j < t < 50t v (finding the values £4 = 1-6 ± 0.1 and 
Ce = 2.±0.1). 

It needs to be emphasized, however, that the limits 
induced by the finiteness of volume and of the inertial 
range extension in both DNS and experimental data do 
not allow for making a definitive statement about the 
behavior in the region r > lOr^. We may ask instead 
if the relative extension of the interval where we see the 
large dip at r ~ 2r, ; and the possible plateau, observed 
for t > 10t i( both in the numerical and experimental 
data (see EXP4 data set), becomes larger or smaller at 
increasing the Reynolds number^ 3 - If the dip region -the 
one presumably affected by vortex filaments- flattens, it 
would give the asymptotically stable scaling properties of 
Lagrangian turbulence. If instead the apparent plateau 
region, at large times, increases in size while the effect of 
high intensity vortex remains limited to time lags around 
(2 -T- 6)t, ; , the plateau region would give the asymptotic 
scaling properties of Lagrangian turbulence. This point 
remains a very important question for the future because, 
as of today, it can not be answered conclusively neither 
by experiments nor by simulations. 



1. Finite volume effects at large time lags 

As noted above, the EXP4 data for C4( T ) develop an 
apparent plateau at a smaller value than the DNS data. 
In this section, we show how the DNS data can be used 
to suggest a possible origin for this mismatch. 

We investigate the behavior of the local slopes for the 
simulations, when the volume of size £ 3 , where parti- 
cles are tracked, is systematically decreased. Essentially 
only trajectories which stay in this sub- volume are con- 
sidered in the analysis, mimicking what happens in the 
experimental measurement volume. We considered vol- 
ume sizes C in the range which goes from the full box 
size B to £>/7, and we average over all the sub-boxes to 
increase the statistical samples. In Fig. El we plot the 
statistics of the trajectory durations for both the experi- 
ment and DNS by varying the measurement volume size. 
For C = £>/4, the modified DNS statistics are essentially 
indistinguishable from the experimental results. It is now 
interesting to look at the LVSF measured from these fi- 
nite length numerical trajectories. This shows that the 




FIG. 6: Comparison of the probability V(t) that a trajec- 
tory lasts a time t vs t/r v for the experiment EXP2 and 
for DNS2 trajectories in different numerical measurement do- 
mains C/B = 1, 1/2, 1/4, 1/6, 1/7. 



method we devised is able to mimic the presence of a 
finite measurament volume as in experiments. 



In Fig. 7(a) 



we show the fourth-order LVSF obtained 
by considering the full length trajectories and the trajec- 
tories living in a sub- volume as explained above. What 



clearly appears from Fig. 7(a) is that the finite length 



of the trajectories lowers the value of the structure func- 
tions for time lags of the order of 20^ < t < 40^. 
Indeed, the finite-length statistics give a signal that is 
always lower than the full averaged quantity: this effect 
may be due to a bias to slow, less energetic particles, 
which have a tendency to linger inside the volume for 
longer times than fast particles, introducing a systematic 
change in the statistics. Note that this is the same trend 
detected when comparing EXP2 and EXP4 in Figs. [31 In 
Fig. |7(b)| we also show the effect of the finite measure- 
ment volume on the local slope for p = 4. By decreasing 
the observation volume, we observe a trend towards a 
shorter and shorter plateau with smaller and smaller val- 
ues. This could be the source of the small offset between 
the plateaux developed by the EXP4 data and the DNS 
data in Fig. [5] 

For the sake of clarity, we should recall that in the DNS 
particles can travel across a cubic fully periodic volume, 
so during their full history they can reenter the volume 
several times. In principle, this may affect the results for 
long time delays. However, since the particle velocity is 
taken at different times we may expect that possible spu- 
rious correlations induced by the periodicity to be very 
small, if not absent. This is indeed confirmed in Fig. |7(b)] 
where we can notice the perfect agreement between data 
obtained by using periodic boundary conditions or lim- 
iting the analysis to subvolumes of size C = B (i.e. not 
retaining the periodicity) and even C = B/2. 
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FIG. 7: (a) The fourth-order structure function £4(1") vs t/t v 
measured from DNS trajectories, for both full length trajec- 
tories (and with periodic boundary conditions) and for trajec- 
tories in smaller measurement volumes C/B = 1/4, 1/7. (b) 
The logarithmic local slope £i(t) measured from DNS trajec- 
tories, for both the full length trajectories (periodic boundary 
conditions) and for trajectories in smaller measurement vol- 
umes C/B — 1,1/2,1/4,1/7. Note the tendency toward a 
less developed plateau, at smaller and smaller values, as the 
measurement volume decreases. 



2. Filtering and measurement error effects at small time 

lags 



FIG. 8: (a) Logarithmic local slope £4(1") for the DNS2 data 
set. The symbol DNS2a,b, and c denote the DNS2 trajec- 
tories modified by noise and filter effects, mimicking what 
was done in the experiments. In particular, DNS2a refers to 
the introduction of noise in the particle position of the order 
of Sx ~ 77/IO and with a Gaussian filter width a ~ t v /3, 
DNS2b to the same filter width but with much larger spa- 
tial noise (Sx ~ fl/A), and DNS2c to the same spatial noise 
but a large filter width a ~ 2r^/3. Note how when the fil- 
ter is not very large and with large spatial errors we have 
strong non-monotonic behavior for the local slopes (DNS2b). 
(b) The effect of filter width on data from EXP2 experiment 
(Rx = 690, small measurement volume). We tested 4 different 
filter widths: a/r v = 1/6, 1/3, 2/3, and 4/3. 



As discussed in Sect. UH results at small time lags can 
be slightly contaminated by several effects both in DNS 
and experiments. DNS data can be biased by resolu- 
tion effects due to interpolation of the Eulerian velocity 
field at the particle position. In experiments uncorre- 
cted experimental noise needs to be filtered to recover 
the trajectories ! 7 ' 11 ^ 6 

To understand the importance of such effects quan- 
titatively, we have modified the numerical Lagrangian 
trajectories in the following way. First, we have intro- 



duced a random noise of the order of 77/IO to the parti- 
cle position, in order to mimic the noise present in the 
experimental particle detection. Second, we have imple- 
mented the same Gaussian filter of variable width used 
to smooth the experimental trajectories x(t). We also 
tested the effect of filtering by processing experimental 
data with filters of different length. 

In Figs. |8(a)| and |8(b)[ we show the local scaling expo- 
nents for £4(7-) as measured from these modified DNS tra- 
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jectories together with the results obtained from the ex- 
periment, for several filter widths. The qualitative trend 
is very similar for both the DNS and the experiment. 
The noise in particle position introduces non-monotonic 
behavior in the local slopes at very small time lags in the 
DNS trajectories. This effect clearly indicates that small 
scale noise may strongly perturb measurements at small 
time lags, but will not have important consequences for 
the behavior on time scales larger than t v . On the other 
hand, the effect of the filter is to increase the smoothness 
at small time lags slightly (notice the shift of local slopes 
curves toward the right for r ~ t v for increasing filter 
widths). A similar trend is observed in the experimental 
data (Fig. |8(b)| . In this case, choosing the filter width 
to be in the range r £ [1/6, 1/3] t,, seems to be optimal, 
minimizing the dependence on the filter width and the 
effects on the relevant time lags. Understanding filter ef- 
fects may be even more important for experiments with 
larger particles, on the order of or comparable with the 
Kolmogorov scale. In those cases, the particle size natu- 
rally introduces a filtering by averaging velocity fluctua- 
tions over its size, i.e., those particles are not faithfully 
following the fluid trajectories. 1 ^ 15 

IV. CONCLUSION AND PERSPECTIVES 

A detailed comparison between state-of-the-art exper- 
imental and numerical data of Lagrangian statistics in 
turbulent flows has been presented. The focus has been 
on single-particle Lagrangian structure functions. Only 
due to the critical comparison of experimental and DNS 
data it is possible to achieve a quantitative understand- 
ing of the velocity scaling properties over the entire range 
of time scales, and for a wide range of Reynolds numbers. 

In particular, the availability of high Reynolds number 
experimental measurements allowed us to assess in a ro- 
bust way the existence of very intense fluctuations, with 
high intermittency in the Lagrangian statistics around 
r £ [2 : 6]t v . For larger time lags r > 10t v , the signature 
of different statistics seems to emerge, with again good 
agreement between DNS and experiment (see Fig. [5]). 
Whether the trend of logarithmic local slopes at large 
times is becoming more and more extended at larger and 
larger Reynolds number is an issue for further research. 

Both experiments and numerics show in the ESS local 
slope of the fourth and sixth order Lagrangian structure 
functions a dip region at around time lags (2-^6)r^ and a 
flattening at r > 10t, ; . As of today, it is unclear whether 
the dip or the flattening region give the asymptotic scal- 
ing properties of Lagrangian turbulence. The question of 



which region will extend as a function of Reynolds num- 
ber can not be resolved at present, and remains open for 
future research. 

It would also be important to probe the possible 
relations between Eulerian and Lagrangian statistics 
as suggested by simple phenomenological multifractal 
modelS i 13 i 23 ' 47 i 48 In these models, the translation be- 
tween Eulerian (single-time) spatial statistics and La- 
grangian statistics is made via the dimensional expression 
of the local eddy turnover time at scale r: r r ~ r / 5 r u. 
This allows predictions for Lagrangian statistics if the 
Eulerian counterpart is known. An interesting applica- 
tion concerns Lagrangian acceleration statistics^ where 
this procedure has given excellent agreement with exper- 
imental measurements. When applied to single-particle 
velocities, multifractal predictions for the LVSF scaling 
exponents are close to the plateau values observed in 
DNS at time lags r > 10t v . It is not at all clear, how- 
ever, if this formalism is able to capture the complex 
behavior of the local scaling exponents close to the dip 
region r G [2 : 6]t v , as depicted in Fig. [5j Indeed, mul- 
tifractal phenomenology, as with all multiplicative ran- 
dom cascade models^ does not contain any signature 
of spatial structures such as vortex filaments. It is pos- 
sible that in the Lagrangian framework a more refined 
matching to the viscous dissipative scaling is needed, as 
was proposed in Ref. [T3 , rephrasing known results for 
Eulerian statistics 4i Even less clear is the relevance for 
Lagrangian turbulence of other phenomenological mod- 
els, based on super-statistics 4 ^, as recently questioned in 
Ref. H3 

The formulation of a stochastic model able to capture 
the whole shape of local scaling properties from the small- 
est to the largest time lag, as depicted in Fig. remains 
an open important theoretical challenge. 
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